Bow Shock Research Articles

Energy Transport and Conversion Within Earth's Supercritical Bow Shock: The Role of Intense Lower‐Hybrid Whistler Waves

AbstractDetailed analysis of a high Mach number quasiperpendicular Earth bow shock crossing by the Magnetospheric Multiscale (MMS) spacecraft fleet reveal that lower‐hybrid (LH) whistler waves generated in the shock foot region transport energy predominately along the shock surface and slightly toward the shock ramp in the shock normal incidence frame, where wave energy accumulates and is dissipated into the plasma. This suggests the LH whistlers play an integral role in energy reconfiguration at high Mach number collisionless shocks with ramifications to plasma heating. The multipoint observations are used to quantify the wave characteristic parameters (via interferometry), Poynting fluxes, and energy conversion rates D, and to assess their scale dependencies and spatial and temporal properties. The whistler associated energy transport and conversion are found to depend on scale and location within the layer. High‐frequency electrostatic waves yield largest values of D. However, the dominant net energy exchange contribution is from the LH whistlers. In the foot spatially temporally coherent net energy exchange from the plasma to whistlers is observed, whereas deeper in the ramp net wave energy dissipation to the plasma is observed exhibiting significant space‐time variability. These results are consistent with the modified two stream instability driven by the relative drift between reflected ions and electrons as the mechanism for wave growth in the foot. Owing to strong electron heating, whistler energy dissipation in the ramp is attributed to Landau damping, which out‐competes the destabilizing effect of the reflected ion and electron drift.

Changes in and Recovery of the Turbulence Properties in the Magnetosheath for Different Solar Wind Streams

Solar wind is known to have different properties depending on its origin at the Sun. In addition to the differences in plasma and magnetic field parameters, these streams differ due to the properties of turbulent fluctuations involved in the flow. The present study addresses the changes in the turbulence properties in the magnetosheath—the transition region in front of the magnetosphere. This study is based on statistics from the simultaneous measurements of magnetic field fluctuations in the solar wind and in the magnetosheath. Both the dayside and flank magnetosheath regions are focused on to detect the evolution of the turbulent fluctuations during their flow around the magnetosphere. Turbulent cascade is shown to save its properties for fast solar wind streams. Conditions favorable for the preservation of the turbulence properties at the bow shock may correspond to the increased geoefficiency of large-scale solar wind structures.

Unusually Low-frequency Whistler-mode Waves and Their Association with High-energy Protons Observed Upstream of Martian Bow Shock

Whistler-mode waves upstream of planetary bow shock are often referred to as “1-Hz waves” due to the center of their observed frequency range being at ∼1 Hz. A series of whistler-mode waves were observed upstream of the Martian bow shock by MAVEN on 2015 August 14, with unusually low frequencies centered at ∼0.4 Hz. These waves were accompanied (though not synchronized) by the significant flux enhancement of high-energy protons up to ∼10 keV. By analyzing the wave dispersion property and the wave–particle interaction condition, we find that the observed whistler-mode waves have the potential of resonating with protons of ∼1 keV with large pitch angles up to nearly perpendicular to the background magnetic field, thereby providing a feasible means of accounting for proton acceleration. Our results indicate the possible origin of energized protons in the Martian environment through the interaction with whistler-mode waves, and their potential relationship with the unique upstream conditions.

Investigating Protostellar Accretion-driven Outflows across the Mass Spectrum: JWST NIRSpec Integral Field Unit 3–5 μm Spectral Mapping of Five Young Protostars

Investigating Protostellar Accretion is a Cycle 1 JWST program using the NIRSpec+MIRI integral field units to obtain 2.9–28 μm spectral cubes of five young protostars with luminosities of 0.2–10,000 L ⊙ in their primary accretion phase. This paper introduces the NIRSpec 2.9–5.3 μm data of the inner 840–9000 au with spatial resolutions from 28 to 300 au. The spectra show rising continuum emission; deep ice absorption; emission from H2, H i, and [Fe ii]; and the CO fundamental series in emission and absorption. Maps of the continuum emission show scattered light cavities for all five protostars. In the cavities, collimated jets are detected in [Fe ii] for the four <320 L ⊙ protostars, two of which are additionally traced in Brα. Knots of [Fe ii] emission are detected toward the most luminous protostar, and knots of [Fe ii] emission with dynamical times of <30 yr are found in the jets of the others. While only one jet is traced in H2, knots of H2 and CO are detected in the jets of four protostars. H2 is seen extending through the cavities, showing that they are filled by warm molecular gas. Bright H2 emission is seen along the walls of a single cavity, while in three cavities narrow shells of H2 emission are found, one of which has an [Fe ii] knot at its apex. These data show cavities containing collimated jets traced in atomic/ionic gas surrounded by warm molecular gas in a wide-angle wind and/or gas accelerated by bow shocks in the jets.

Quantification of magnetosphere–ionosphere coupling timescales using mutual information: response of terrestrial radio emissions and ionospheric–magnetospheric currents

Abstract. Auroral kilometric radiation (AKR) is a terrestrial radio emission excited by the same accelerated electrons which excite auroral emissions. Although it is well correlated with auroral and geomagnetic activity, the coupling timescales between AKR and different magnetospheric or ionospheric regions have yet to be determined. Estimation of these coupling timescales is non-trivial as a result of complex, non-linear processes which rarely occur in isolation. In this study, the mutual information between AKR intensity and different geomagnetic indices is used to assess the correlation between variables. Indices are shifted to different temporal lags relative to AKR intensity, and the lag at which the variables have the most shared information is found. This lag is interpreted as the coupling timescale. The AKR source region receives the effects of a shared driver before the auroral ionosphere. Conversely, the polar ionosphere reacts to a shared driver before the AKR source region. Bow shock interplanetary magnetic field BZ is excited about 1 h before AKR enhancements. This work provides quantitatively determined temporal context to the coupling timelines at Earth. The results suggest that there is a sequence of excitation following the onset of a shared driver: first, the polar ionosphere feels the effects, followed by the AKR source region and then the auroral ionosphere.

Statistical analysis of high-frequency whistler waves at Earth's bow shock: Further support for stochastic shock drift acceleration

We statistically investigate high-frequency whistler waves (with frequencies higher than ∼10% of the local electron cyclotron frequency) at Earth's bow shock using magnetospheric multi-scale (MMS) spacecraft observations. We focus specifically on the wave power within the shock transition layer, where we expect electron acceleration via stochastic shock drift acceleration (SSDA) to occur associated with efficient pitch-angle scattering by whistler waves. We find that the wave power is positively correlated with both the Alfvén Mach number in the normal incidence frame MA and in the de Hoffmann–Teller frame MA/cos θBn. The empirical relation with MA/cos θBn is compared with the theory of SSDA that predicts a threshold wave power proportional to (MA/cos θBn)−2. The result suggests that the wave power exceeds the theoretical threshold for MA/cos θBn≳30–60, beyond which efficient electron acceleration is expected. This aligns very well with previous statistical analysis of electron acceleration at Earth's bow shock [Oka et al., Geophys. Res. Lett. 33, 5–6 (2006)]. Therefore, we consider that this study provides further support for SSDA as the mechanism of electron acceleration at Earth's bow shock. At higher-Mach-number astrophysical shocks, SSDA will be able to inject electrons into the diffusive shock acceleration process for subsequent acceleration to cosmic-ray energies.

Cross‐Comparison of Observations With the Predictions of Global Hybrid Simulations for Multiple IMF Discontinuities Impacting the Bow Shock and Magnetosheath

AbstractWe use the three‐dimensional (3‐D) global hybrid code ANGIE3D to simulate the interaction of four solar wind tangential discontinuities (TDs) observed by ARTEMIS P1 from 0740 UT to 0800 UT on 28 December 2019 with the bow shock, magnetosheath, and magnetosphere. We demonstrate how the four discontinuities produce foreshock transients, a magnetosheath cavity‐like structure, and a brief magnetopause crossing observed by THEMIS and MMS spacecraft from 0800 UT to 0830 UT. THEMIS D observed entries into foreshock transients exhibiting low density, low magnetic field strength, and high temperature cores bounded by compressional regions with high densities and high magnetic field strengths. The MMS spacecraft observed cavities with strongly depressed magnetic field strengths and highly deflected velocity in the magnetosheath downstream from the foreshock. Dawnside THEMIS A magnetosheath observations indicate a brief magnetosphere entry exhibiting enhanced magnetic field strength, low density, and decreased and deflected velocity (sunward flow). The solar wind inputs into the 3‐D hybrid simulations resemble those seen by ARTEMIS. We simulate the interaction of four oblique TDs with properties similar to those in the observation. We place virtual spacecraft at the locations where observations were made. The hybrid simulations predict similar characteristics of the foreshock transients, a magnetosheath cavity, and a magnetopause crossing with characteristics similar to those observed by the multi‐spacecraft observations. The detailed and successful comparison of the interaction involving multiple TDs will be presented.

Simulating Wind-Blown Nebulae from Single and Binary Massive Stars

Winds from massive stars expand supersonically into their surroundings, creating dynamic and fascinating nebulae that can give us insight into physical processes in interstellar plasma, and into the evolutionary history of the stars. Around single stars, parsec-scale bubbles such as bow shocks and ring nebulae are formed, whereas in colliding-wind binary (CWB) systems the high wind density produces intense time- and space-dependent emission across the electromagnetic spectrum from radio to gamma-rays. This contribution summarizes some recent results from 3D MHD modelling of bow shocks around runaway stars such as ζ Oph, and of the wind-collision zone of the CWB systems WR140 and WR21a. A resolution study of 3D simulations of bow shocks shows that X-ray emission from the shocked wind is time-variable and that converged results can be obtained once the Kelvin-Helmholtz instability at the contact discontinuity is resolved. Simulations of the CWB system WR140 show that inverse-Compton cooling of the shocked plasma can trigger runaway cooling when the orbit is near periastron, producing strong compression and dynamical instabilities. This sharply reduces the hard-X-ray emission around periastron, in agreement with observations. Scaling tests of the simulation software pion are also presented for a model of the CWB system WR21a run on up to 8192 cores using the HPC system Karolina.

Tracking the Subsolar Bow Shock and Magnetopause: Applying the Magnetosheath Velocity Gradient Method

AbstractBoth the bow shock and magnetopause move in response to varying solar wind and magnetospheric conditions. Tracking their locations can provide important clues to the state of the solar wind‐magnetosphere interaction, but is difficult with single spacecraft observations. This paper employs multipoint THEMIS observations of velocity gradients in the subsolar magnetosheath to remotely sense boundary locations on a continual basis for various solar wind conditions. We present three cases: (a) continuous northward IMF and no magnetopause motion; (b) southward IMF and no magnetopause motion with evidence of nightside activity; and (c) southward IMF and pressure increase with inward motion. When observing spacecraft are located near the Sun‐Earth line, inferred boundary locations agree well with the predictions of the BATS‐R‐US global magnetohydrodynamic model, confirming the utility of both the new method and the models. Results show that boundaries often lie nearly at rest with amplitudes less than 0.5 RE. They provide evidence indicating that nightside reconnection and a strong sunward convection in the outer magnetosphere can counteract the magnetopause erosion expected when a southward interplanetary magnetic field (IMF) initiates reconnection on the dayside magnetopause.

Experimental study on control of transverse jet mixing by arrayed plasma energy deposition

The efficient and prompt mixing of fuel is crucial in the operation of scramjet engines. This paper presents the findings from wind tunnel experiments that examined the influence of plasma energy deposition on transverse jets at a Mach number of 6.13. The study took into account various inlet flow total pressures and momentum flux ratios between the jet and the main flow. Utilizing a database containing time-resolved intensities from instantaneous schlieren images, we perform turbulence analysis employing various techniques such as the root mean square, fast Fourier transform, proper orthogonal decomposition, and the two-point correlation method. Specifically, we aim to compare and analyze the pulsation characteristics and spatial self-organization of the jet flow field, both with and without energy deposition control. The findings reveal that intermittent “hot bubbles” created by plasma energy deposition interact with the bow shock induced by the jet, resulting in the formation of an array of large-scale vortices. These vortices emerge as the dominant structures within the jet, effectively amplifying its pulsations. At low inlet flow pressures, energy deposition primarily disrupts the jet, causing large-scale vortices to propagate primarily within the jet plume region. However, at high inlet flow pressures, the impact of energy deposition extends to both the jet and the turbulent boundary layer, encompassing their respective disturbance ranges. Increasing the inlet flow pressure constraints the evolution of large-scale vortices, thus limiting the efficacy of energy deposition in governing the mixing process.

FAUST

Context. The origin of the chemical diversity observed around low-mass protostars probably resides in the earliest history of these systems. Aims. We aim to investigate the impact of protostellar feedback on the chemistry and grain growth in the circumstellar medium of multiple stellar systems. Methods. In the context of the ALMA Large Program FAUST, we present high-resolution (50 au) observations of CH3OH, H2CO, and SiO and continuum emission at 1.3 mm and 3 mm towards the Corona Australis star cluster. Results. Methanol emission reveals an arc-like structure at ∼1800 au from the protostellar system IRS7B along the direction perpendicular to the major axis of the disc. The arc is located at the edge of two elongated continuum structures that define a cone emerging from IRS7B. The region inside the cone is probed by H2CO, while the eastern wall of the arc shows bright emission in SiO, a typical shock tracer. Taking into account the association with a previously detected radio jet imaged with JVLA at 6 cm, the molecular arc reveals for the first time a bow shock driven by IRS7B and a two-sided dust cavity opened by the mass-loss process. For each cavity wall, we derive an average H2 column density of ∼7 × 1021 cm−2, a mass of ∼9 × 10−3 M⊙, and a lower limit on the dust spectral index of 1.4. Conclusions. These observations provide the first evidence of a shock and a conical dust cavity opened by the jet driven by IRS7B, with important implications for the chemical enrichment and grain growth in the envelope of Solar System analogues.

Ion Dynamics Across a Low Mach Number Bow Shock

AbstractA thorough understanding of collisionless shocks requires knowledge of how different ion species are accelerated across the shock. We investigate a bow shock crossing using the Magnetospheric Multiscale spacecraft after a coronal mass ejection crossed Earth, which led to solar wind consisting of protons, alpha particles, and singly charged helium ions. The three species are resolved upstream of the shock. The low Mach number of the bow shock enabled the ions to be partly distinguished downstream of the shock due to the relatively low ion heating. Some of the protons are specularly reflected and produce quasi‐periodic fine structures in the velocity distribution functions downstream of the shock. Heavier ions are shown to transit the shock without reflection. However, the gyromotion of the heavier ions partially obscures the fine structure of proton distributions. Additionally, the calculated proton moments are unreliable when the different ion species are not distinguished by the particle detector. The need for high time‐resolution mass‐resolving ion detectors when investigating collisionless shocks is discussed.

Accuracy of Global Geospace Simulations: Influence of Solar Wind Monitor Location and Solar Wind Driving

AbstractSome space weather models, such as the Space Weather Modeling Framework (SWMF) used in this study, use solar wind propagated from the first Lagrange point (L1) to the bow shock nose (BSN) to forecast geomagnetic storms. The SWMF is a highly coupled framework of space weather models that include multiple facets of the Geospace environment, such as the magnetosphere and ionosphere. The propagated solar wind measurements are used as a boundary condition for SWMF. The solar wind propagation method is a timeshift based on the calculated phase front normal (PFN) which leads to some uncertainties. For example, the propagated solar wind could have evolved during this timeshift. We use a data set of 123 geomagnetic storms between 2010 and 2019 run by the SWMF Geospace configuration to analyze the impact solar wind propagation and solar wind driving has on the geomagnetic indices. We look at the probability distributions of errors in SYM‐H, cross polar cap potential (CPCP), and auroral electrojet indices AL and AU. Through studying the median errors (MdE), standard deviations and standardized regression coefficients, we find that the errors depend on the propagation parameters. Among the results, we show that the accuracy of the simulated SYM‐H depends on the spacecraft distance from the Sun‐Earth line. We also quantify the dependence of the standard deviation in SYM‐H errors on the PFN and solar wind pressure. These statistics provide an insight into how the propagation method affects the final product of the simulation, which are the geomagnetic indices.

Aerodynamic Design and Performance Analysis of Mars Ascent Vehicles

The mid-lift-to-drag ratio vehicle possesses larger dimensions and enhanced aerodynamic characteristics compared to the conventional sphere-cone reentry vehicle, suggesting its suitability for crewed missions to Mars. This research employs Computational Fluid Dynamics (CFD) simulations to investigate the influence of varying geometries on the aerodynamic properties of mid-lift-to-drag ratio vehicles. The aerodynamic analysis of four distinct configurations demonstrates that a bow shock wave forms on the windward side, leading to a reduction in flow velocity and subsequent increases in pressure and temperature after the shock wave. The simulation results for the four configurations in both laminar and turbulent flow conditions, specifically using the Shear-Stress Transport (SST) model, reveal that turbulence results in higher aerodynamic heating compared to laminar flow. Moreover, the aerodynamic heating is greater at a 40° angle of attack than at 30°. In this study, the transition phenomena of four geometric configurations are examined using the k-ω-γ transition model. The findings indicate a pronounced transition to turbulence on the windward side of the Ellipsled 2.00-0.25 and Hammerhead-Nominal configurations under conditions of high Reynolds numbers, transitioning from laminar to turbulent flow. The investigation into the aerodynamic characteristics of mid-lift-to-drag ratio vehicles holds substantial importance for ensuring the safety and success of Martian reentry missions, highlighting the critical nature of aerodynamic behavior in the context of space exploration.

Shocks, clouds, and atomic outflows in active galactic nuclei hosting relativistic jets

Context. A number of observations have revealed atomic and/or molecular lines in active galaxies hosting jets and outflows. Line widths indicate outward motions of hundreds to a few thousands of kilometers per second. They appear to be associated with the presence of radio emission in Gigahert-peaked spectrum (GPS) and/or compact steep spectrum (CSS) sources, with linear sizes of ≤10 kpc. Numerical simulations have shown that the bow shocks triggered by relativistic jets in their host galaxies drive ionization and turbulence in the interstellar medium (ISM). However, the presence of atomic lines requires rapid recombination of ionized gas, which appears challenging to explain from the physical conditions revealed thus far based on numerical simulations of powerful jets. Aims. The aim of this paper is to provide a global framework to explain the presence of lines in terms of jet and shock evolution and to fix the parameter space where the atomic and molecular outflows might occur. Methods. This parameter space is inspired by numerical simulations and basic analytical models of jet evolution as a background. Results. Our results show that a plausible general explanation involves momentum transfer and heating to the interstellar medium gas by jet triggered shocks within the inner kiloparsecs. The presence of post-shock atomic gas is possible in the case of shocks interacting with dense clouds that remain relatively stable after the shock passage. Conclusions. According to our results, current numerical simulations cannot reproduce the physical conditions to explain the presence of atomic and molecular outflows in young radio sources. However, I show that these outflows might occur in low-power jets at all scales and I predict a trend towards powerful jets showing lines at CSS scales, when clouds have cooled to recombination temperatures.

Quantified Plasma Heating and Energy Dissipation in the Earth’s Quasi-perpendicular Bow Shock

The high-speed solar wind decelerates as it travels through the Earth’s bow shock, during which a portion of the incident flow energy is converted into the plasma thermal energy. However, the energy partition of plasma heating among different species, as well as the energy dissipation mechanism, remains elusive. In this paper, we quantify the plasma heating and the energy dissipation by calculating the pressure–strain interactions in 33 quasi-perpendicular bow shocks observed by the Magnetospheric Multiscale spacecraft. Our results show that the Joule dissipation measured by J·E′ and the quasi-viscous dissipation quantified by PiD are distinct in the examined shocks. We also reveal that ions gain more energy than electrons, and the compressive effect is more important than the incompressible channel in plasma heating. PiD contributed by the gyrotropic pressure tensor is consistently positive for electrons; however, there is no discernible distinction between the gyrotropic and nongyrotropic contribution to PiD for ions.

The nature of the X-ray filaments around bow shock pulsar wind nebulae

Context. We propose that the X-ray filaments emerging from selected bow shock pulsar wind nebulae are due to a charge-separated outflow of electrons and/or positrons escaping the nebula and propagating along the local Galactic magnetic field. Aims. The X-ray brightness, length, and thickness of filaments are all accounted for if a nonresonant streaming instability is excited. Methods. This is possible if particles are released in the interstellar medium as a collimated beam, as would be expected in a reconnection region between the nebular and interstellar magnetic fields. Results. We successfully test this idea on the Guitar Nebula filament and discuss other cases. Conclusions. These filaments provide the best diagnostics available for particle escape from evolved pulsar wind nebulae, a process essential to assessing the contribution of these sources to cosmic ray positrons. The same phenomenology might govern the occurrence of TeV halos and their importance for cosmic ray transport.

Investigation of multi-scale flow structures and combustion characteristics in a cavity-enhanced circular scramjet

This study delves into the investigation of multi-scale flow structures and combustion characteristics in a cavity-enhanced circular scramjet operating under a Mach-4.5 freestream flow condition with a total temperature of 2200 K and a total pressure of 164 kPa. To achieve this, a hybrid large-eddy/Reynolds-averaged Navier-Stokes method, coupled with a pressure-corrected Flamelet/Progress Variable model, is employed and validated. The combustion flow field shows abundant small- and medium-scale flow structures induced by compressibility and combustion heat release. Vorticity alteration is found to be mainly caused by these factors, with the volumetric expansion term and baroclinic term serving as the primary driving terms. The whole combustion process is characterized by zoning and multi-scale features in spatial, temporal, heat release and species distribution, arising from differences in fuel-air mixing state and chemical reaction rates in the flow field. The injection-formed bow shock and the cavity ramp-formed compression wave interact intensively with the fuel plume, enriching the turbulent structures and strongly coupling with the chemical reaction process. The cavity not only promotes fuel-air mixing and combustion, but also facilitate the rapid transition from the laminar-like flame to thickened turbulent combustion. Statistical analysis reveals that the combustion heat release is predominantly governed by the scramjet mode and the diffusion flames, and further demonstrates the multi-scale flow and chemical reaction characteristics. Most of the combustion flow field is observed to be governed by fast chemistry, primarily in the corrugated flamelet regime, while some slow chemistry zones are also identified, lying in the distributed reaction zone regime.

Stellar wind bubbles of OB stars as Galactic cosmic-ray re-accelerators

Abstract Cosmic rays are highly energetic messengers propagating in magnetized plasma, which are, possibly but not exclusively, accelerated at astrophysical shocks. Amongst the variety of astrophysical objects presenting shocks, the huge circumstellar stellar wind bubbles forming around very massive stars, are potential non-thermal emitters. We present the 1D magneto-hydrodynamical simulation of the evolving magnetized surroundings of a single, OB-type main-sequence $60\, \rm M_{\odot }$ star, which is post-processed to calculate the re-acceleration of pre-existing non-thermal particles of the Galactic cosmic ray background. It is found that the forward shock of such circumstellar bubble can, during the early phase ($1\, \rm Myr$) of its expansion, act as a substantial re-accelerator of pre-existing interstellar cosmic rays. This results in an increasing excess emission flux by a factor of 5, the hadronic component producing γ-rays by π0 decay being more important than those by synchrotron and inverse Compton radiation mechanisms. We propose that this effect is at work in the circumstellar environments of massive stars in general and we conjecture that other nebulae such as the stellar wind bow shocks of runaway massive stars also act as Galactic cosmic-ray re-accelerators. Particularly, this study supports the interpretation of the enhanced hadronic emission flux measured from the surroundings of κ Ori as originating from the acceleration of pre-existing particles at the forward shock of its wind bubble.

A nozzle design for mitigating particle slowing in the bow shock region during micro-cold spray of 8 YSZ films

During micro-cold spray deposition, also referred to as the aerosol deposition method or vacuum kinetic spraying, performed with conventional nozzles, the particle impact velocities decrease drastically with particle size for fine particles <500 nm in diameter due to slowing in the stagnation region. A new design for a nozzle that contains pressure relief channels is proposed that allows the pressure in the stagnation region downstream of the bow shock to be reduced. This reduced stagnation pressure results in less particle slowing compared to conventional nozzle geometries, particularly for smaller or less dense particles. The effects of the channel geometry on the particle impact velocity are systematically investigated by independently varying the channel parameters. Calculations show that the impact velocities for 100 and 200 nm yttria stabilized zirconia particles is increased by 111% and 31%, respectively, for a selected pressure relief channel nozzle when compared to a comparable conventional nozzle. Although impact velocities are increased, a tradeoff exists with this nozzle design in that the particle focusing in the nozzle is decreased and some of the particles may be removed from the aerosol by the channels. Experiments using a nitrogen as a carrier gas at ∼40 kPa upstream pressure show that, despite the loss of larger particles into the relief channels, the deposition efficiency is improved by 300% when depositing films from fine 8 YSZ powder with the pressure relief nozzle compared to a conventional nozzle.